Enhanced sensitivity to time- variation of m p /m e in the inversion spectrum of ammonia 
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We calculate the sensitivity of the inversion spectrum of ammonia to possible time-variation 
of the ratio of the proton mass to the electron mass, yi = m p /m e . For the inversion transition 
(A ~ 1.25 cm -1 ) the relative frequency shift is significantly enhanced: 8lu/uj — — 4.46<5/i//i. This 
enhancement allows one to increase sensitivity to the time- variation of \i using NH3 spectra for high 
redshift objects. We use published data on microwave spectra of the object B0218+357 to place 
the limit Sfi/fi = (0.6 ± 1.9) x 10~ 6 at redshift 2 = 0.6847; this limit is several times better than 
the limits obtained by different methods and may be significantly improved. Assuming linear time 



dependence we obtain p,/fi = ( — 1 ± 3) x 10 1(1 
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INTRODUCTION 



The possible time-variation of the fundamental con- 
stants has been discussed for a long time. The interest 
in this discussion has grown considerably after the re- 
cent discovery of the acceleration of the universe. The 
latter is usually regarded as evidence for the existence 
of dark energy. Cosmological evolution of dark energy 
may cause variations in fundamental constants, such as 
the fine-structure constant a and the proton to electron 
mass ratio, fi = m p /m e . The electron mass is one of 
the parameters of the Standard Model, it is proportional 
to the vacuum expectation value of the Higgs field (the 
weak scale). The proton mass is proportional to another 
fundamental parameter, the quantum chromodynamics 
scale Aqqd {m p ~ 3Aqcd)- The proportionality coef- 
ficients cancel out in the relative variation. Therefore, 
we are speaking about the relative variation of a very 
important dimensionless fundamental parameter of the 
Standard Model, the ratio of the strong to weak scale, 
defined as 6(AQ C D/m e ) / '{h.QCD/m e ) = Sfi/fi. 

It is known that /i defines the scales of electronic, vi- 
brational, and rotational intervals in molecular spectra, 
E c \ ■ E vih : E rot ~ 1 : pT 1 ! 2 : jjT 1 . Similarly, the ratio of 
electronic and hyperfine intervals in atoms and molecules 
also depend on /1, E c \ : £Jhfs ~ 1 : « ffpM ) where g p is the 
proton g-factor. These scalings are used to look for the 
time-variation of /z by comparing electronic, vibrational, 
rotational, and hyperfine spectra of atoms and molecules 
P, [i| . In the most recent astrophysical studies Q a non- 
zero effect was reported for two quasars at 3.5 a level: 



Sfi/fi = (20 ±6) x 10" 



(1) 



at a time scale of approximately 12 Gyr. Assuming lin- 
ear variation with time this result translates into \xj \x = 
(— 17±5)xl0~ 16 yr~ 1 . A different method, compari- 
son of the hyperfine transition in atomic hydrogen with 



optical transitions in ions, was used in Refs. 0, 0|. This 
method allows one to study variation of the parameter 
x = a 2 g p / \i. Analysis of 9 quasar spectra with redshifts 
0.23 < z < 2.35 gave 



Sx/x = (6.3 ±9.9) x 10~ 6 , 
x/x = (-6 ±12) x 10~ 16 



yr" 



(2) 
(3) 



which is consistent with zero variation of /1. In Refs. 6, 
0,Q the 18 cm A-doublet lines in OH molecule were stud- 
ied from objects at the redshifts z « 0.247, z ks 0.6847, 
and z w 0.765 and no time-variation of the parameter 
g p {c? [i) v was seen, where v < 2. 

In this paper we use the enhanced sensitivity of the 
inversion spectrum of ammonia to variation of fi to place 
a new limit on the time- variation of /z at the cosmologi- 
cal timescale. The NH3 molecule has a pyramidal shape 
and the inversion frequency depends on the exponentially 
small tunneling of three hydrogens through the potential 
barrier Because of that it is very sensitive to any 
changes of the parameters of the system, particularly to 
the reduced mass for this vibrational mode (loj . 

We use high-resolution ammonia spectra for gravita- 
tional lens B0218±357, published by Henkel et al. [H|. 
The redshifts for ammonia lines are compared to the red- 
shifts for the rotational lines of other molecules measured 
in Refs. [Ill Il2t 1 1 31 ] - The ammonia lines have an order 
of magnitude stronger dependence on /1 than the usual 
vibrational lines; this enhancement allows us to place the 
best limit on the variation of /x. 



INVERSION SPECTRUM OF NH 3 

The inversion spectrum of NH3 has been studied for a 
very long time [9j and is considered a classical example 
of the tunneling phenomenon. The inversion vibrational 
mode is described by a double well potential with first 



2 



two vibrational levels lying below the barrier. Because of 
the tunneling, these two levels are split in inversion dou- 
blets. The lower doublet corresponds to the wavelength 
A ~ 1.25 cm and is used in ammonia masers. Molec- 
ular rotation leads to the centrifugal distortion of the 
potential curve. Because of that the inversion splitting 
depends on the rotational angular momentum J and its 
projection on the molecular symmetry axis K: 



r(J,K)=> 



Cl [j{J+l)-K 2 ]+c 2 K 2 , (4) 



where we omitted terms with higher powers of J and K. 
Numerically, ujf nv ps 23.787 GHz, a ps 151.3 MHz, and 
c 2 ps 59.7 MHz. 

In addition to the rotational structure ((4|) the inver- 
sion spectrum includes much smaller hyperfine structure. 
For the main nitrogen isotope 1 N, the hyperfine struc- 
ture is dominated by the electric quadrupole interaction 
(~ 1 MHz) [3]. Because of the dipole selection rule 
AK = the levels with J — K are metastable and in 
laboratory experiments the width of the corresponding 
inversion lines is usually determined by collisional broad- 
ening. In astrophysics, the hyperfine structure for spectra 
with high rcdshifts is not resolved and we will not discuss 
it here. 

For our purposes it is important to know how the 
parameters in (j4|) depend on fundamental constants. 
One can measure only dimcnsionlcss ratios of frequencies 
which do not depend on the units used. It is convenient 
to consider all parameters in atomic units. The energy 
unit Hartree is Eh = m e e A /h 2 = e 2 /dB, where as is the 
Bohr radius (Eh— 2 Ry=219475 cm -1 ). In these units 
all electron energies (E e /Ejj) and electrostatic poten- 
tials (U(r)/Ejj) have no dependence on the fundamental 
constants (here we neglect small relativistic corrections 
which give a weak a dependence), the vibrational inter- 
vals ~ yT 1 ! 2 and the rotational intervals ~ A* -1 - The 
inversion frequency uif DV / Eh and constants ci^/Eh are 
also functions of fi only (see below). Note that the coef- 
ficients Ci depend on /i through the reduced mass of the 
inversion mode and because they are inversely propor- 
tional to the molecular moments of inertia. That implies 
a different scaling of u)f av and c; with fi. The magnetic 
hyperfine structure of NH3 is due to the interaction of 
nuclear magnetic moments and proportional to a 2 g 2 fi~ 2 . 

We see that different frequencies in the inversion spec- 
trum scale differently with \i and a. In principle, this 
allows one to study time-variation of /j, and a by com- 
paring different lines of the inversion spectrum. On the 
other hand, it may be preferable to use independent ref- 
erences (see below). 



INVERSION HAMILTONIAN 

The inversion spectrum Q can be approximately de- 
scribed by the following Hamiltonian: 



H'uiv — 2M1 ®x 



2M 

+ T^[J(J + 1)-K 2 ] 



(5) 



K 



where x is the distance from N to the H-plane, I\, I2 
are moments of inertia perpendicular and parallel to the 
molecular axis correspondingly and Mi is the reduced 
mass for the inversion mode. If we assume that the length 
d of the N — H bond does not change during inversion, 
then Mi — 2.54m p and 



\m v d 2 [l 



h(x) 
l2(x) ~ 3m p d 2 [l 



-0.2{x/d) 2 ] 
(x/d) 2 ] . 



(6) 
(7) 



The dependence of Ii t 2 on x generates a correction to the 
potential energy of the form C(J, K)x 2 / '/j,. This changes 
the vibrational frequency and the effective height of the 
potential barrier, therefore changing the inversion fre- 
quency Winv given by Eq. Q . 

Following [15j we can write the potential U(x) in (O 
in the following form: 



U(x) — ^kx 2 + &exp (— ca 



(8) 



Fitting vibrational frequencies for NH3 and ND3 gives 
k ps 0.7598 a.u., b ps 0.05684 a.u., and c ps 1.3696 a.u. 
Numerical integration of the Schrodinger equation with 
potential ([8]) gives the following result: 



SuJu 



-4.46 



5(i 



(9) 



It is instructive to reproduce this result from an analyt- 
ical calculation. In the semiclassical approximation the 
inversion frequency is estimated as (l6j : 



(10a) 



x /2Mi(U(x) - E)dx , (10b) 



^inv = — exp (-5") 

LU V ( 1 

" 6XP (~h 



where uj v is the vibrational frequency of the inversion 
mode, S is the action in units of h, x = ±<z are classical 
turning points for the energy E. For the lowest vibra- 
tional state E = Umin + \lo v . Using the experimental 
values 0^=950 cm -1 and cc>i nv =0.8 cm" 1 , we get S~5.9 . 

Expression (|10b|) allows one to calculate 
the dependence of uj® nv on the mass ratio 
fi. Let us present S in the following form: 
S = An 1 ' 2 /" o y/{U(x) - E)/E H d(x/a B ), where A 
is a numerical constant. We see that the dependence of 
w hiv 011 I 1 appears from the factor [i 1 / 2 in S and from 
the vibrational frequency u> v and E — U m [ n = \>-o v which 



3 



are proportional to [i^ 1 ! 2 . Below we assume that all 
energies are measured in atomic units and omit the 
atomic energy unit E^j. Then we obtain 

1 dS dSdE\ , 
2^ + d^ + dE^)> (llb) 

where we took into account that dS/da = because the 
integrand in (|10b[) turns to zero at x = ±a. 

It is easy to see that dS/dn = S/2[i. The value of 
the third term in Eq. (|llb[) depends on the form of the 
potential barrier: 

dS q S 



dE 



4C/„ 



E 1 



where for the square barrier q = 1, and for the triangular 
barrier q = 3. For a more realistic barrier shape q k 2. 
Using paramctrization (|8|) to determine Z7 ma x we get: 

W 1 + 5+ g yT J) =-4.4^. (13) 



2a* V ^ ^max - £-7 A* 

We see that the inversion frequency of NH3 is an order 
of magnitude more sensitive to the change of fj, than typ- 
ical vibrational frequencies. The reason for that is clear 
from Eq. (TT3"|) : it is the large value of the action S for the 
tunneling process. 

Let us also find the dependence of the constants on 
(i in Eq. (g]). According to Eqs. © - ([7]) both constants 
must have the same dependence on \x. Below we focus 
on the constant C2, which is linked to the last term in 
Hamiltonian It follows from Eq. ([7]) that this term 
generates a correction to the potential: 



SU(x) = 



K 2 



3m p cfi 



(14) 



This correction does not change the height of the barrier, 
but changes the energy E = C/ m i n + \w v in (|10b|) by rais- 
ing the potential minimum and increasing the vibrational 
frequency: 



U n 



K 2 



3m p c? 4 



'■<)■■ 



1 



K 2 



(15) 
(16) 



With the help of Eq. (|T2|) with q = 2 we can find the 
constant C2: 



C2 



3m p d 4 k 



1 



4=s 



E 



(17) 



We can differentiate Eq. ifTTj) to estimate how C2 depends 
on fx. This leads to: 802/02 = — 5.0<5aVa*j while the nu- 
merical solution with Hamiltonian ([5]) gives: 



Sci :2 _ 1 Sfl 

= —5.1 — . 

C1.2 A* 



(18) 



It is clear that NH3 is not the only molecule with en- 
hanced sensitivity to variation of [i. Similar enhancement 
should take place for all tunneling transitions in molecu- 
lar spectra. For example, the inversion frequency for ND 3 
molecule is 15 times smaller than for NH3 and Eq. QlOap 
leads to S « 8.4, compared to S w 5.9 for NH3. Accord- 
ing to Eq. (|T3|) that leads to a slightly higher sensitivity 
of the inversion frequency to \x [1 01 ] : 



ND, 



8C2 
C2 



-5.7^ 

1" 



-6.2 ^ 



(19) 



REDSHIFTS FOR MOLECULAR LINES IN THE 
MICROWAVE SPECTRA OF B0218+357 

In the previous section we saw that the inversion fre- 
quency ujf nv and the rotational intervals w mv (Ji,Ki) — 
Wi nv (J2, K2) have different dependencies on the constant 
Ai. In principle, that allows one to study time- variation of 
A* by comparing different intervals in the inversion spec- 
trum of ammonia. For example, if we compare the rota- 
tional interval to the inversion frequency, then Eqs. ([9]) 
and l[T8|) give: 



6{[cJi nv (Ji,Ki) 

UJ'mv 

(J 2 ,K 2 )]/u? nv } 

[Winv(Jl, Ki) - L0- mv (J 2 , K 2 )}/w° nv 



.0.6^. (20) 
A* 



The relative effects are substantially larger if we compare 
the inversion transitions with the transitions between the 
quadrupole and magnetic hyperfine components. How- 
ever, in practice this method will not work because of the 
smallncss of the hyperfine structure compared to typical 
line widths in astrophysics. 

It is more promising to compare the inversion spectrum 
of NH3 with rotational spectra of other molecules, where 



^rot 
Wrot 



Sfj, 
At 



(21) 



In astrophysics any frequency shift is related to a corre- 
sponding redshift: 



5u> 



1 



(22) 



According to Eqs. j9|) and (f2Tj) . for a given astrophysical 
object with z = zq variation of [i will lead to a change of 
the redshifts of all rotational lines 8z Tai — (1 + zo)dfi/fi 
and corresponding shifts of all inversion lines of ammo- 
nia 5z{ nv = 4.46 (1 + zq) Sfj,/ fi. Therefore, comparing the 
redshift for NH3 with the redshifts for rotational lines we 
can find Sfj,/ fj,: 



5 ± = 0.289 Ziav 



•2-rot 



1 + zq 



(23) 
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TABLE I: Redshifts for molecular rotational lines, ammonia 

inversion lines, and hydrogen hyperfine line in the spectrum 
of B0218+357. 

Rotational lines 

CO J = 1 -> 2 red-shifted 0.68470 [12] 

blue-shifted 0.68463 [12] 

CO, HCO+, HCN average 0.68466(1) [17] 

Inversion lines of NH3 

NH3 (J,K) = (1,1) red-shifted 0.684679(3) [11] 

blue-shifted 0.684649(15) [11] 

= (2, 2) red-shifted 0.684677(3) [11] 

blue-shifted 0.684650(17) [11] 

= (3, 3) red-shifted 0.684673(3) [11] 

blue-shifted 0.684627(33) [11] 
average red-shifted 0.684676(3) 

average blue-shifted 0.684647(11) 

H A = 21 cm average 0.68466(4) [18] 



In Table U we list the redshifts for microwave lines in 
the spectrum of the object B0218+357. Three inver- 
sion lines (J, K) = (1, 1), (2, 2), and (3, 3) are reported in 
Ref. [ll[ • Each of them consists of a narrow red-shifted 
and a wide blue-shifted component. The splitting be- 
tween the red-shifted and blue-shifted components, which 
is about 5 km/s, is ascribed to the complicated struc- 



ture of the molecular cloud ll| . Using average red- 
shifts of these inversion components (0.684676(3) and 
0.684647(11)) from Table U we can calculate the aver- 
age deviation of the inversion redshift in respect to the 
average molecular redshift (0.68466(1)): 



a unweighted 



= (0.2 ±0.9) x 10~ 5 , (24) 
Az wei g hted = ( 6 ± 9 ) x 10 - 5 ( 25 ) 



Eq. 



gives the following estimate for variation of /i: 



6 ± = 10- 6 x 



0.3 ± 1.6 (unweighted), 
1.1 ±1.5 (weighted). 



(26) 



As a final result we present a conservative limit with 
larger error bars to cover the total interval between the 
minimal and maximal values for both estimates: 



— = (0.6 ±1.9) x 10~ 6 . 



(27) 



We can also compare averaged redshift for ammonia 
with that of hydrogen to get a restriction on the variation 



of the parameter y = a g p fi : 



try 
V 



^hfs 



1 



Z0 



= (1 ± 17) x 10" 



(28) 



The estimates (|26H28|) can be further improved by 
dedicated analysis of the molecular spectra published in 
As it was mentioned in [121 ] . the major- 



Refs. 11, 12 



com pon ents. The same applies to the hydrogenic 21 cm 
[l9(. Instead of taking an average, as we have done in 



line 



all red-shifted and all blue-shifted components 
should be analyzed independently. That may allow one 
to reduce the error bars significantly. 
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